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In our previous work we have found a lower bound for the multipartite uncertainty product of 
the position and momentum observables over all separable states. In this work we are trying to 
minimize this uncertainty product over a broader class of states to find the fundamental limits 
imposed by nature on the observable quantites. We show that it is necessary to consider pure states 
only and find the infimum of the uncertainty product over a special class of pure states (states 
with spherically symmetric wave functions). It is shown that this infimum is not attained. We 
also explicitly construct a parametrized family of states that approaches the infimum by varying 
the parameter. Since the constructed states beat the lower bound for separable states, they are 
entangled. We thus show that there is a gap that separates the values of a simple measurable 
quantity for separable states from entangled ones and we also try to find the size of this gap. 
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I. INTRODUCTION 

The famous Heisenberg uncertainty relation states that 
the standard deviation of position and momentum ob- 
servables of any quantum states cannot be simultane- 
ously small or, more precisely, that (in appropriate units) 



a x a p ^ 



(1) 



This inequality is tight and becomes equality for the vac- 
uum state whose wave function in the position represen- 
tation is 



1 



(2) 



In pQ we have shown that the inequalities for differ- 
ent degrees of freedom can be "multiplied" side-by-side 
to produce inequalities valid for all completely separable 
states. For example, the inequalities 



and 
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(3) 



are valid for all bipartite separable states. These inequal- 
ities for separable states are also tight since for the two- 
dimensional vacuum they become equality. We have also 
shown that these inequalities can be violated. A question 
of fundamental interest is what is infimum of the prod- 
ucts on left hand side of the inequalities ^ over the set 
of all quantum states? 

For any physical quantum state these products are 
strictly positive, but this does not mean that the infimum 
is also positive. The problem of determining this infimum 
is surprisingly more complicated than in a single partite 
case. In this work we do not completely solve this prob- 
lem, we restrict it to the class of states with real spher- 
ically symmetric wave functions and find the infimum 
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over this class of states. This infimum is not attained, 
so we construct family of states parametrized by a real 
parameter < £ < 1 such that when £ — > 1 the uncer- 
tainty product approaches the infimum. More precisely, 
we construct a family of states such that cr xp a px < 1/4 
for all < £ < 1 and a xp a px — » 1/8 when £ — > 1. We also 
generalize our construction to a larger number of parties. 
As the number of parties growth, the problem becomes 
more and more complicated, so we obtain explicit results 
only for four-partite and six-partite cases and outline the 
general approach. 

The paper is organized as follows. In section [TT] we 
set up the environment in which we work in this paper 
and prove that to minimize the uncertainty product un- 
der consideration pure states are sufficient. In section III 
we derive an analytical expression for the state obtained 
numerically in our previous work and introduce the tech- 
nique that is used to obtain the main result of this paper. 
In section [TV| we obtain our main result, i.e. we derive 
family of bipartite states that minimize the uncertainty 
product in the special class of pure states with real spher- 
ically symmetric wave functions. The derivation is not 
difficult but lengthy, so it is divided into several steps to 
make it easier to follow. In section |V| we analyze some 
properties of these states. In section [Vl| we generalize our 
construction to a more general case of a larger number 
of parties. The conclusion is in section VII. Proofs of the 
auxiliary results are given in the appendices at the end 
of the paper. We make an intensive use of the book [2J, 
which we refer to simply as GR. It is rather a large book 
with a lot of results, so we give the page and formula 
numbers to easily locate the results we refer to. 



II. MULTIPARTITE UNCERTAINTY PRODUCT 

The general form of the uncertainty relation states 
that for arbitrary observables A and A' and for all quan- 
tum states the inequality ua^A 1 ^ 1/2| ([^4, A'])\ is valid, 

where a a = \ (A 2 ) — (A) 2 is the standard deviation of 
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the observable A. It has been shown in Q] that all com- 
pletely separable TV-partite states satisfy the inequality 

cta^.a^a',...^ > ^\([A U [A N ,A' N })\ (4) 

for arbitrary observables Ai and A\ acting on the ith 
part, i = 1, . . . , N. If we take Ai — Xi, A[ — pi for all i, 
we then get the inequality 
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This inequality is tight — it is easy to verify that for the 
iV-partite vacuum state with the wave-function 



ip (xi,. . . ,X N ) 



1 
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(6) 



the left-hand side of the inequality ^ is exactly 2 
The natural question is — how strong can the inequality 
([5} be violated? 

First of all note that by analogy with the inequality 
5|) the following 2 JV ~ 1 inequalities are also valid: 
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(7) 



where the two sequences of x and p complement each 
other — if one of them has x in some position then the 
other has p in the same position. There are 2 N sequences 
of x and p of length N, but the number of different in- 
equalities is only half of this number, i.e. 2 N ~ 1 . If we 
can find a state that violates one inequality of this kind 
then we can find states that violate any such inequality. 
To do it we need the phase-shifting operator [3] 



U v = e~ X{fm = exp 



2 { dx? 



(8) 



The action of this operator on wave functions is simply 
the fractional Fourier transform [4] 
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As one can easily see, if <p — ir/2 then 



(9) 



ip{x) dx. 



J> /2 [^(x)](p) 



^(x)e~ ipx dx = 4>{p) (10) 



R 



is the momentum representation of the state \ijf) with the 
wave function ip(x) in the coordinate representation. Let 
us take a pure bipartite (N = 2) state with the wave func- 
tion ip(x,y) and apply partial fractional Fourier trans- 
form J r n /2 to this wave function with respect to the sec- 
ond argument. Denote this new wave function ip. Then 
for the new state we have the relations 

&xp — 0~ XX i ®px — @pp ■ ( 1 1) 



This means that if the original state violates the inequal- 
ity Q for N = 2 then the new state violates the inequal- 
ity 



> 1/4. 



(12) 



The similar conclusion is also valid in the general mul- 
tipartite case — from a state that violates one of the 
inequalities Q we can construct 2 N ~ 1 — 1 states that 
violate the other 2 N ~ 1 — 1 inequalities of this form. 

If we want to minimize the inequalities Q then we can 
restrict our attention to pure states only due to the 



Theorem 1. If the inequality 



VAB&A'B' ^ S > 



(13) 



is valid for all bipartite pure quantum states of some 
quantum system, where operators A, A' act on one de- 
gree of freedom and B , B' act on the other one, then 
this inequality is also valid for all states (i.e. including 
mixed states). The similar statement can be extended to 
multipartite case. 

The statement of this theorem is intuitively clear be- 
cause by mixturing quantum states we can only increase 
the dispersion. Equivalently it can be formulated as fol- 
lows: 



inf a abo- A' B' 

all 



= inf (JAB^A'B 1 
pure 



(14) 



where all means all bipartite quantum states and pure 
means bipartite pure states only. So, if we want to mini- 
mize the inequalities ([7|) then we should focus our efforts 
on pure states. It is a very useful result since a wave func- 
tion is a much simpler object than a density operator (it 
is not necessary to care about positivity). 

Proof. The proof follows the idea of Q]. We have the 
inequality 



2 

a AB 



^ 2 °~A'B' ^ 2\o-abo~a'B' ^ 2A<5, 



(15) 



which is valid for all bipartite pure states and for all 
non-negative A by assumption. Now let us take a general 
mixed state g. It can be written as a mixture of pure 
states 



= ^2pk\ipk)(^k\ 



(16) 



where pu are non-negative numbers such that ^2 k Pk = 1 
and \4>k) are some bipartite (pure) states. Due to the 
concavity of the dispersion, we have 



°~ab{q) + >?°~A>B>{e) 

fe 

> 2X5j2Pk = 2X6. 



(17) 
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From Lemma 3 of pQ we conclude that 
46 2 ^4<t ab (q)<t 2 a , b ,(q), 



(18) 



which is equivalent to the inequality (|13j) for the mixed 
state g. □ 

As a measure of violation of the inequalities Q we 
take the ratio of the right-hand to the left-hand side. In 
P] a state violating the inequality ( 12 ) has been con- 
structed. The violation of this inequality for that state is 
w 1.2192. In this work we try to do better and find the 
minimal value of the product a xp a px for a special class of 
states, thus constructing highly entangled (with respect 
to this inequality) states and generalize our construction 
for larger numbers of parties. 



III. SIMPLE STATE 



c n = £ n c%, n J? 1 and then express the quadratic form 
Q as a function of Co and Ci only 



Q = -co Cl +cj(2-OYl ( 2 + "K 



- C oc 1 + (2-0( » ^^)cl 



i-e 



(22) 



On the other hand, from the relation for the transformed 
coefficients we immediately obtain the expression for the 
original coefficients c n 



c n = Cl, n ^ 1. 

n 



(23) 



The first two coefficients cq and c\ are not independent, 
from the normalization condition of the state ( 19 1 we 
obtain the following relation between them: 



We start our discussion of states that violate the in- 



equality ( 12 ) with the derivation of an analytically ex- 



pression for the coefficients of the state 



E 

71=0 
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(19) 



that has been used in [T] and get an exact analytical 
expression for the violation value of 1.2192 that has been 
previously obtained numerically. It has been shown that 
the difference a xp a px — 1/4 for this state is given by the 
following quadratic form of the coefficients c n : 



Q 



(2n(2n + l)t& - (n + l)(2n + 1) ). (20) 



This quadratic form has been minimized by numerically 
computing the minimal eigenvalue A m i n ~ —0.04495 of 
the truncated matrix corresponding to this form. The 
components of the normalized eigenvector corresponding 
to the minimal eigenvalue arc the coefficients of the state 



(19) and for this state we have a xp a px = (1/4) + A 



0.20505. We now derive an analytical expression for these 
coefficients. The method presented here for this simpler 
case will be useful later when we present technique to 
solve more complicated and more general problems. 

Factoring out the common terms under the summation 
sign in the form ( 20 ) , we arrive to a simpler expression 



Q = -coci + J2 [ 2 + - ) S «( 2£ ™ - 5 «+l)' ( 21 ) 



where we use the coefficient transformation c„ = nc n , 
n ^ 1. We do not transform the coefficient Cq and the 
coefficient c\ is left unchanged: ci = c%. Let us assume 
that the transformed coefficients satisfy the simple rela- 
tion c„+i = n 1, where |£| < 1 is a real parameter 
to be determined, i.e. that they form a geometrical pro- 
gression. From this relation we can easily derive that 



+°° f2n-2 t : ( t2\ 

4 + clYj—=c 2 + I ^§dcl = l, (24) 



n=l 
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where Li s (z) is the polylogarithm special function defined 
by the infinite power series as [5] 



(25) 



The equation ( 24 ) describes an ellipse and this ellipse can 



be parametrized as follows: 



Cq = COS if, Cl 



: sin if. 



(26) 



Substituting these expressions into the equation (22), we 



obtain the following expression for our quadratic form: 

Q = -Ci(0 sin(2^) + C 2 (0(1 - cos(2^)), (27) 
where the coefficients Ci(£) and C2(£) are given by 
1 £ 



<M0 



C 2 (0 = 



2yLW) : 

(2-0 



m(i-a\ e 



(28) 
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The minimal value of the expression (27) when £ is fixed 
and ip varies is 

Qo(0 - C 2 (0 - ^(O + CfCfl, (29) 

and the value of the angle ip for which this minimal value 
is attained is determined from the equation 



tan(2<ys) 



Ci(0 



(30) 



The plot of the function ( 29 ) is shown in Fig. [T] One 
can numerically minimize this function (for example, 
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FIG. 1: The function Qo(£)> defined by the equation (291 



with Mathematica) and get that the minimal value is 
re -0.04495 for £ min » 0.318674. This result is in perfect 
agreement with the previous result obtained numerically 
by computing the eigenvalues of the quadratic form ( 20 1 . 



To obtain the coefficients c n , we take the parameter £ m j n 
and find the corresponding angle ip nl - m from the equation 
(30 1. Then, we compute the coefficients cq and c\ ac- 



cording to the equation ( 26 1 . The rest of the coefficients 
are obtained from the equation ( 23 ) . Though we have 
not strictly proved that —0.04495 is the minimal eigen- 



value of the quadratic form (20), we have analytically 



constructed a state on which this value is attained and 
shown that this result agrees with the numerical compu- 
tations. 



IV. 



CONSTRUCTION OF SPHERICALLY 
SYMMETRIC STATES 



The state ( 19 ) has a very special form, so one can ask 
whether it possible to stronger violate the inequality ( 12 ) 
with a more general state. Here, we construct such a fam- 
ily of states that violate this inequality. The construction 
is lengthy and thus divided into several steps. 

First step. At this step we analyze the product o~ xp o~ px 
and transform it into a more manageable form. The vari- 



ances <7™ and a 2 read as 



J 2 

xp 



vlp = (%lPb) - (XaPb) 2 , <J 2 VX = (f a xl) - (p a X b ) 2 . (31) 

Their product would be easier to deal with if the averages 
(x a Pb) and (p a Xb) were zero. This is the case when, for 
example, the wave function ip(x, y) of the state under- 
study is real. In fact, for a real wave function we have 

ft dip 

(x a Pb) = -i J J ^ip{x,y)x — (x,y)dxdy = 0, (32) 

since for any fixed x the integral over y can be taken 
explicitly 



^ dy = l^ j2{x > +M) - ^ = °' (33) 



and is equal to zero due to the normalization property of 
ip. In the same way we obtain the equality (p a Xb) = 0. 



We have that for a real wave function we can write 

2 2 1-2 -2\ / -2 »2\ 

VxpVpx = ( x aPb){Pa x b) 

1 (34) 

1 1 1 -2 -2 i "2 -2\2 1-2-2 -2 -2\2\ v ' 
= ^{( X aPb +Pa X b) - ( X aPb ~ Pa X b) ) 

and conclude that the product a xp a px can be bounded 
from above by the square root of the first term on the 
right-hand side of this expression, i.e. the following in- 
equality is valid: 

<7x P <J px < ^( x IpI +pI x I) = (%)■ (35) 

So, at this step, we have reduced our minimization prob- 
lem to the study of the operator Z on the set of real wave 
functions. 

Second step. In the standard position representation, 
the operator Z is given by the following differential op- 
erator: 



V l/~2~2 i -2 -2\ 1 ( 2 

Z = T^aPb +Pa x b) = [ X Qy2 



dx 2 



(36) 



The average value of this operator on a real wave function 
ip(x,y) can be computed as 



(Z) 



R 2 



V> I x 2 



dy 2 



' dxdy. 



I) 



dx 2 



(37) 



Integrating by parts, this expression can be transformed 
as follows: 



d 2 V> , dip 

WjT2 dx = V7T 
R ax z ox 



+oo 



R 



dip 
dx 



dx. 



(38) 



For a normalizable wave function the first term on the 
right-hand side is zero and, transforming the other term 
under the integral in Eq. (37), we can conclude that for 
the average (Z) we have 



(Z) = 




fa) 



dx dy. (39) 



Looking at this expression, it seems natural to write it 
in polar coordinates. The derivatives with respect to the 
Cartesian coordinates x and y can be expressed in terms 
of the derivatives with respect to polar coordinates r and 
6 as follows: 



d_ 

dx 
0_ 
dy 



= sm 



dr 
dr 



sin 8 



cos 9 



ld_ 
rd~6' 
ld_ 
rd~6' 



(40) 



If we substitute these derivatives into the integral on 
the right-hand side of Eq. (39 1, the resulting expression 



will not look any simpler, so we will consider only real 
spherically symmetric wave function ip(x,y), i.e. func- 
tions that do not depend on the angle 9 when written 
in polar coordinates as ^(r, 9) = ip(r cos9,r sin 9). Then 
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#(r,0) = *(r), so d^/dO = and Eq. ([39]) takes a sim- 
pler form 



(Z) 



1 



^ ; / r 2 sin 2 (26») f ^- ) ,W 



o Jo 
o 



'f I ,1, 

dr 



(41) 



The normalization of the real wave function ip{x,y) can 
be written as follows: 

/■+0O />27T 

^ 2 (x,y)dxdy = / / ^ 2 (r)rdrd9 

R 2 JO JO Oq*) 

/> + 00 V ' 

= 2tt / ty 2 (r)rdr = 1. 
Jo 

Here we note that the second term on right-hand side of 

, d 2 



the equation ( 34 1 is 



-2-2 -2-2 
X aPb Pa X b 



dy 2 



y 



d 2 

dx 2 ' 



(43) 



and for our spherically symmetric wave function ip(x,y) 
we have 



R- 



dy 2 



dx 2 



dx dy 



ijj x 



2 / d i> 



R 2 



dy 




(44) 



<ix c??; = 0, 



which can be easily obtained with the help of Eqs. (401 



This means that in the case of a real spherically symmet- 
ric wave function we simply have 



(Z). 



(45) 



In this step, our original problem has been further re- 
duced to the problem of minimizing the integral (41 1 pro- 



vided that the function ^(r) satisfies the normalization 
condition (42). 



Third step. To simplify the integral (41) and the nor- 



malization condition ( 42 ) let us introduce the function 
f(r) via the relation 



*(r) 



1 



f(r 2 )- 



(46) 



The normalization condition ( 42 1 in terms of this func- 
tion reads as 

2?r / ^ 2 (r)rdr = 2tt -f(r 2 )rdr 

Jo Jo * (47) 



f{r 2 )dr 2 = 



t(r)dr = l. 



We see that the function /(r) is normalized in the ordi- 
nary sense. Now let us write the last integral of Eq. (41 ) 



in terms of the function f(r). For the derivative we have 



dr 



-^rf(r 2 ), 



(48) 



and, substituting this expression into Eq. ( 41 ) , we get 

2 



(Z) = 







dr J 



dr= r°f 2 (r z )dr 



\J r^f 2 (r 2 )dr 2 = \^ r 2 f 2 (r)dr. 



(49) 



In this step, we have formulated our problem as follows: 
find the minimum of the integral 



(Z) = 



1 



r 2 f' 2 (r)dr 



provided that the function /(r) is normalized 



-t-oo 



f{r)dr=l. 



(50) 



(51) 



In other words, we need to minimize the functional (50) 



over the set of normalized functions. Intuitively it seems 
to be clear that the minimum is not attained because 



function for which the integral ( 50 ) is small should have 
the shape of a peak located near the origin r = 0. Out- 
side the peak the function should be small so that it is 
nearly constant and thus the derivative is almost zero, 
and the more narrow the peak, the smaller the region 
where f'(r) is nonzero. On the other hand, the function 
/ must be normalized, and the more narrow the peak 
the higher it must be, and then one can expect the the 
derivative inside the peak is large. The range of r where 
f'(r) is nonzero is small but the value of the derivative is 
large, so that the function must have a special shape for 
the combination of these two competing features to give 
the smallest value of the integral (50). So, we expect to 



find a parametrized family of functions which represent 
infinitesimally narrow and infinitely high peaks of some 
special form when the parameter varies. 

Fourth step. Let us show that the value of the right- 
hand side of Eq. (50) cannot be smaller than 1/8. To 
prove this, we first compute the integral 



J rf(r)f(r)dr = -rf(r)X°° 



+ 00 



f 2 (r)dr 



(52) 



If we apply the Cauchy-Schwarz inequality 

^f(r)9(r)dr\ < J* f 2 (r) dr j\ 2 {r) dr, (53) 

which is valid for all real integrable functions /(r) and 
g(r) on an interval [a, b], — oo ^ a < b ^ +oo, to the 
functions f(r) and g(r) = rf'{r) on the interval [0, +oo), 
we get 



f 2 (r)dr 



r 2 f' 2 (r)dr, 



(54) 



G 



and the normalization condition (51 1 gives us the desired 
inequality 



1 f + °° I 



which is valid for n ^ (for n = assume that i_i(r) = 
0). From this equality, for a general function (60) one 
can easily obtain the relation 



(55) 



+00 



Note that, in fact, we have just proved the inequality 



Ik/T > i||/|| 2 , 



(56) 



where we use the standard notation for the norm and 
scalar product of two square-integrable functions /, g <E 
L(0,+oo) 



rf'(r) = -^^(-nd n -i+d n + (n+l)d n+1 )l n (r). (63) 

Using the orthonormality of {l n ( r )} one can get 

/ r 2 f 2 (r) dr = - V(rf„ + (n + l)d n+1 - nd n -if 
J o 4 n=0 



r+oo 

(/,<?) = / f(r)g(r)dr, 
Jo 

r+oo 

ll/f = (/,/)= / f(r)dr. 
Jo 



(57) 



= \ + ^(^0^1,^2, ■ ■ ■), 



where R = R(do, d%, d 2 , . . .) is the quadratic form 



(64) 



This notation and the inequality ([56} will be used for the R = £ ( n ( n + x )^n - (n + l)(n+ 2)d n d n+2 ) . (65) 



multipartite generalizations of the results obtained here. 

The Cauchy-Schwarz inequality ( |53[ ) becomes equality 
if and only if the functions f(r) and g(r) are linearly de- 
pendent on the interval [a, b\. In our case this condition 
reads as rf'(r) + A/(r) = for all r ^ 0. The general 
solution of this equation is f(r) — Cr~ x . One can easily 
see that such a function is not normalizable on the in- 
terval [0, +oo). This means that there is no normalized 
functions f(r) that minimize the inequality (55). The 
problem now is to construct such functions f(r) that ap- 
proach to the lower bound 1/8 as close as possible. As a 
hint in this construction we look for functions for which 
there is a linear combination of rf'(r) and f(r) that is 
small in some sense. 

Fifth step. Using the Laguerre functions 



n=0 



From the inequality (55) we can conclude that R ^ —1/2 



provided that the equality (61) is satisfied. Going back 



to Eqs. (45) and (49), we have 
2 r+°° 



Cxp&px — \Z) — 



1 1 



r{r)dr = - + -R. (66) 



l n {r)=L n (r)e- r ' 2 , 



0,1,.. 



(58) 



We thus reduced our minimization problem to the mini- 
mization problem of the quadratic form ( |65| ). 

Sixth step. In this step we construct a parametrized 
family of states for which R tends to —1/2, and thus 
VxpVpx tends to 1/8, but this minimum is not attained. 
Computing the eigenvectors corresponding to the mini- 
mal eigenvalue of the truncated matrices of the quadratic 
R one can observe that the components with odd indices 
are zero, so let us set dm+\ — 0, d 2n — c n , n ^ 0. Then 



which are complete and orthonormal on the interval 
[0, +oo) [B], where L n (r) are Laguerre polynomials [7] 



the form (65) becomes 



p r d, n 



(59) 



we can write any function f(r), normalizable on [0, +oo), 
as a linear combination of the Laguerre functions 



R = J2(M^ + 1)<£ - 2(n + l)(2n + l)c„c„ +1 ). (67) 



This is very similar to the quadratic form Q defined by 



Eq. ( 20 ) , except for the additional factor 2 in the second 



/(r) = 5>»Z n (r). 



(60) 



The normalization of f(r) and orthonormality of {l n {f)} 
give us the condition 



term under the summation sign. One can try the same 
approach we used for the form Q, but in this case it does 
not give us the optimal state, so another idea is needed. 
From Eq. (63) we can obtain 



f'(r) + \f{r) = \ ^(nd n _! - (n + l)d n+1 )l n (r 



(61) 



n=0 



From the relations for Laguerre polynomials [S] one can 
derive the following equality: 

<(r) = ~(nZ n _!(r) + l n (r) - (n + (62) 



* +oo 

= - 5^((2n + l)c„ - (2n + 2)c n+1 )l 2n+1 (r). 

n=0 

(68) 

If the relation c n+ i/c n = (2n + l)/(2n + 2) were valid 
for all n ^ 0, then we would have linear dependence of 
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rf'(r) and f(r) but, as we know, it is impossible for a 
normalized function f(r). And, in fact, one can easily 
check that the sequence c„, defined by this ratio, cannot 
be normalized. But we can try to define the coefficients 
via 



Cn+l 



2n + 1 
2n + 2 



(69) 



where < £ < 1. From this relation one can derive an 
explicit expression for the coefficient c„ 

{2n\ /A" (2n- 1)!' „ , , 

for n ^ 1. The normalization reads as 

/ +00 

h + E 



^ ' (2n-l)!! \ 2 2 r \ 2 2 

v ~(2^irJ C Jco--TOcb-l, (71) 



where if (£) is the complete elliptic integral of the first 
kind E 



jr/2 



(10 



2n 



1 *e 

16 



I 1+ E 



- - (2n-l)!! \ 2 2 „ N 



We thus obtain the coefficient c 



Co 



2if(0' 



(72) 



(73) 



From Eq. ( 68 ) we get 



rf'{r) + \f{r) = LA 



+ 00 



^(2n+l)c n i 2n+1 (r), (74) 



n=0 



and using orthonormality of {l n ( r )} again we derive 



+00/ 1 n2 , 1-^2 ±5° 

(r/ / (r) + -/(r)J * = 



n=0 



2g(Q - (1 - g 2 )if (Q 
4(1+£ 2 M£) ' 



(75) 



where is the complete elliptic integral of the second 
kind [TO] 



tt/2 



l-£ 2 sin 2 9d9 



+ OO 



(2n-l)!!V £ 
(2n)!! 



c2;j 



(76) 



2n-l 



Since lim^! = 1 and lim^! if (£) — +00, we can 
conclude that the integral on the left-hand side of the 



equality (751 tends to zero when £—>•!. One can say 



that the linear combination rf'(r) + (l/2)/(r) becomes 
smaller and tends to zero when £ — > 1. This means that 
our choice of coefficients ( 69 1 is a good candidate for the 
minimizing function. 



Now we can write the function f(r) defined by Eq. ( 60 ) 
as (remember that e?2n+i — and d 2n = c n) 



/(r) = 5> n L 2 „(r) e -''/ 2 



71 = 



+ OO 



i'-E 2 i ,( )>• ! - 

\ n—l 



n 



(77) 



where cq is given by Eq. ( 73 ). To simplify the sum in this 



expression, note that the Laguerre polynomial L n (r) can 
be represented as follows [XT] : 



Ln(r) = e' 



+00 , n 



J (2Vrt)e~ t dt, 



(78) 



where Jq(z) is the Bessel function of the first kind [T2] 



+00 



JbW = E(- 1 ) B ss 



n=0 



2 2 «(n!) 2 ' 



(79) 



Substituting the integral representation |78| ) into Eq. ( 77 ) 
we get 



(80) 



f(r) = c e r / 2 / t) Jo(2Vrt)e-* d£ 



f" /2ra-l\ t 2 ™ 



where S(£,t) is the sum 



+00 



s(^) = i + E^bi 



ra=l 



2n / (2n)! 



(81) 



Note that we can simplify this sum as follows: 

s(t t) = i + T 1 ^-^^^ 

' ^ 2 2 "- 1 n!(n- 1)! (2n)! 



+ OO 



(82) 



where io(.z) = Jo(iz) is the modified Bessel function of 
the first kind. Finally, we arrive at the following expres- 
sion: 

/ — — — r+00 

/(r) = v 21m &T ' 2 Jo I °^ J ^ 2 ^ e ~ t dt 

In Appendix X we prove that the exchange summation 
and integration in Eq. (77) is legal. In this step, we have 



constructed a family of the functions (83), parametrized 



by the real parameter < £ < 1, that are good can- 
didates for the functions that approach the equality in 



Eq. (55) as £ ^ 1. 
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FIG. 2: The quantity R denned by Eq. (871 



Seventh step. Now let us compute the quadratic form 



(67). Substituting the coefficients (70) into the expres- 



sion (67), we get R = (R\ — i?2)co, where Ri and i?2 are 
defined via 



Ri 



£g (2n + l)H (2n-l)!! 2n 
^ (2n)!! (2n-2)!P ' 



+ OO 



n=0 



(2n + l)H 
(2n)!! 



(84) 



The first sum can be simplified as follows: 



ri=0 



(2«; 



ii 



2)1 



^ 2 (i+a^)-(i-e 2 )^(o 
(i-e 2 ) 2 

The second sum can be computed analogously 



(85) 



R 2 = £ 



E 

n=0 



(2n- 



1) 



7T V <2£ 
= 2|(2£(0 

7T 



(2n-l)!! 
(2n)!! 



2)1 



(86) 



We finally obtain 

R = 



(87) 



i+e (i + o 2 ^(0' 

As Fig. [2] illustrates, J? < for all < £ < 1. From the 
expression ( |87| ) we also see that R — > —1/2 when £ — >■ 1. 
In this step we have proved that with the choice of the 



coefficients ( 69 ) the form ( 67 1 tends to its minimal value 



-1/2, but does not attain it. 
Eighth step. Having proved that the expression (jS3h 



gives us the desired family of functions, we now simplify 
this expression and transform it into a proper integral. 
To do it, we need the following standard theorem that 
can be found, for example, in [20] . 



Theorem 2. Let F(t,9) be a continuous function of two 
arguments on the set [£o,+oo) x [floi^i]- If the integral 

p+ca 

F(9) = / F{t,6)dt (88) 

J t 

converges uniformly on the interval [6q, 9±\, then the func- 
tion F(9) is integrable on this interval and the following 
equality is valid: 



f+oo 

F(9)d9= I d9 I F(t,9)dt 

9 Q J t 

dt I F(t,9)d9. 



(89) 



In other words, under these conditions one can change 
the order of integration. 

Note that the Bessel function Iq(z) can be represented 
as follows [15] : 



1 



(90) 



The Theorem 2 can be used to substitute this expression 



into the equation ( 83 ) and exchange the order of integra- 
. In 

f(r) 



tion. In fact, we have 
1 



,r/2 



+ 00 



^mf'-k H™"' (91) 



where 



F(t,9) = e-^ 1+ ^ cos ^ t J (2VH). 



(92) 



To show that the integral (88) converges uniformly on 



the interval [0, n] for any fixed < £ < 1 note that 

F(t, 9) = e -(i-H»»fl)V£t e -(i-VBt J (2Vrt). (93) 

The function e^ 1- ^* Jo(2\/rt) is integrable on [0, +oo) 
due to the equality [T4"] 



- at J (bVi) dt = 



,-b 2 /(4a) 



(94) 



valid for a > 0. For a = this function is not integrable. 
We have that the integral 



^- ( - 1 -^ t J {2Vrl)dt 



(95) 



converges uniformly with respect to 9, since it does not 
depend on 9, and the function g-^+^^Vit [ s uniformly 



bounded by 1. We conclude that the integral (88) con 



verges uniformly and exchanging the order of integration 
is legal. We then obtain the following expression for the 
function f(r) = f^(r): 



1 



f ( (r) = -^==e T ' 2 I d9 I F(t,9)dt 



d9. 



1— sJT cos 9 r 
I r CX P l - 1 + 7? cos 9 2 



(96) 



y/2irK(£) Jo 1 + VI cos 9 
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This integral (in fact, even more general one) can be 
found in tables |15j . and we have 



2^(0(1-0 Vi -i 



' 1-5 2 



(97) 



Let us show that the this function has finite right deriva- 
tive at the origin. We differentiate the expression (97) 



with respect to r and take the limit of this derivative for 
r — >• 0. This limit is then the right derivative /+(0). We 
have 



/+(0) 



1 



8K (0 7(WF' 



(98) 



and this expression is finite for all < £ < 1. 

In this step we finally obtain the expression for the 
family of wave functions: 



r™ exp (- 

l 



1 



V 



1 

exp (-j 



-VI cose x 2 +y 2 
1+Vfcose 2 



1+| x'+y* 
£ 2 



^2^(0(1-0 



Io 



1 + V^cosd 



d9 



(99) 



(x 2 + y 2 ) 



From the finiteness of the right derivative / + (0) we can 
conclude that this wave function is smooth at the origin, 
since the derivative df(r 2 )/dr = 2rf'(r 2 ) is equal to zero 
for r = 0. The function ^(r), defined by Eq. (46), is 
shown in Fig. [3j In fact, it has meaning only for ^ r < 
+oo, but here it is illustrated for the whole real line. Note 
that the shape of these functions is as we have expected. 




FIG. 3: The wave function \I/(r), defined by Eqs. (46 1 and 
(97 1, for a few different values of £. 



Using the integral representation of the wave function 
( 99 1 , we can compute the product a xp a px according to 



the first equality of Eq. ( 34 ) . Exchanging the orders of 



integration as we did before, we can arrive to the follow- 
ing expression for this product: 



1 



8ttK(0 



*■ (1 — £ cos 2 6){1 — £cos 2 9') 
(1-£cos6>cos6»') 3 



d9d9'. 
(100) 



This expression can be shown to agree with Eqs. (66) 
and (87). In this way we get directly from the definition 



the desired property of this wave function. The bipar- 
tite state with the wave function (99) we denote as 
For all < £ < 1 these states violate the inequality ( [12] ) 
and when £ — > 1, we have a xp a px — > 1/8, so the viola- 
tion tends to 2. In next sections, we study some simple 
properties of the states constructed and then extend this 
approach to a more general multipartite case. 



V. PROPERTIES 

In this section we compute the scalar product of the 
states with different values of the parameters £ and 
find the decomposition of these states in the Fock basis. 



A. Scalar product of states with different 
parameters 

On can easily verify that the scalar product | ) of 
the states |^) and can be computed as the scalar 

product of their corresponding functions /j (r) and fe> (r) 



+ OG 



ft(r)f 6 ,(r)dr. 



(101) 



For the latter we have 

-+oo 2 r+oo 



de 

i) J0 



exp 



1— \/£ cos 9 r 
l + \/Ccos0 2 



l + V^cosS' 2 



(l + V^cos6»)(l + v/fcosfl' 



de 1 . 



(102) 



In this expression the integrations over the angles 9 and 
0' are performed first and then the integration over r. 
Note that 



1 fl - V£cos9 1 - cos(9' 

2 \l + Vfcosfl + 1 + Vfcosfl' 

1 - % /cJFcos6»cos6'' 



(1 + COS0)(1 + y/^COs9' 



(103) 



(i + VD(i + VS0 



> 0. 



From this we can wc can conclude that the order of the 
integrations can be changed and the integration over r 
can be performed first and we get the equality 



+oo 1 

d9d9' 



(104) 



1 



/{£ cos 9 cos 9' 
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After the integration over 8 we get 
+00 1 

d8' 



(105) 



V 1 - W cos2 



Making substitution z = cos 6', the last integral is trans- 
formed to the following one: 



dz 



dz 



-1 ^{l-^'z 2 ){l-z 2 ) ^0 V(l-^ 2 )(l-« 2 ) 

_ ( 106 ) 

which, in turn, is transformed to 2K (V?? 7 ) by the sub- 
stitution z = sin 8. We then obtain the desired scalar 
product 



k(VW) 



(107) 



Note that for £ = £' we get the normalization condition 
= 1, as it must be. 



B. Fock states representation 

Now we derive the coefficients c nm , n,m ^ 0, of the 
state l^f) in the Fock basis. To do this, we need to find 
the scalar product of the state with the Fock state 
\nm) , c nm = (nm|^) or, in other words, we need to 
compute the integral 



^ n {x)ijj m {y)i}){x,y)dxdy, (108) 



R 2 



where 



4>n{x) = . * H n (x)e- x2 ' 2 (109) 
y/y/Tr2 n n\ 

is the wave function of the Fock state \n). To compute 
this integral, we introduce the generating function 

n ,n„,m 

F(x, y, U, v) = lpn(x)lp m (y) 



>/nlm\ 



1 / x 2 + y 2 + u 2 + v 2 r-. , 

exp h V2 m + yv) 

A V 2 



(110) 



This generating function is useful because it is much eas- 
ier to compute the integral 



I(u,v) 



R 2 



F(x,y,u 1 v)ip{x,y)dxdy, (111) 



and then expand it in u and v to find the coefficients c„ 
We have 



I(u,v) 



e 2 



+00 



rdr 



2tt 



(112) 



/(r 2 )e~ r2 / 2 gV / 2r(ii cos 6»+u sine) 



The inner integral (over 8) can be easily taken by noting 
that u cos 8 + vsind = \Ju 2 + v 2 cos(# + U)O ) and that 
due to periodicity of cosines function the shift 8 U , V plays 
no role in the integration over the period: 



I(u,v) = 2e~ 

r-\-oo 







/(r 2 )e- r / 2 / (r v /2(u2 + u 2))rdr = (n3) 



+ OO 



f(r)e- r / 2 I {^2r{u 2 +v 2 )) dr. 



We can substitute the expression (961 for the function 



f(r) and get a repeated integral for I(u,v). One can 
easily check, as it has been done before, that the order 
of integration can be changed and at the end we get the 
following result: 



I(u, v) 



1 



(114) 



The integral can be taken explicitly and expressed in 
terms of the Bessel function I , but the expression (114) 



is more convenient for our purpose. The exponent under 
the integral is the product of two exponents, one contain- 
ing u and the other containing v. We can expand them 
in u and v respectively, multiply and integrate over 8. As 
one can see, we have to integrate powers of the cosines 
function. For these powers we have (see Eq. ( 127) below) 



cos™ 8d8 = - , 

2" \n/2 



(115) 



if n is even and the integral of odd powers are zero. From 
this, one can obtain 



+00 



n + m\ u v 



1 1 m. I 



(116) 

where the sum runs over all n and m with n + m even. 



On the other hand, from Eqs. ([108]) , ([110]) and ([111]) we 

have 



I(u, v) 



^ u n v m 

/ 4 Cnm / , , ' 

n.m— 



(117) 



Comparing this with the expansion (116), we obtain the 
coefficients c„m 



n-\-m 



, (118) 



where either n — An' , m = Am' or n — An' + 2, in = Am' + 
2, all other coefficients are zero. As an additional check, 
in Appendix [B] we derive directly from the expression 
( |118[ ) that the coefficients c nm satisfy the normalization 
condition. 
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VI. GENERALIZATION 

In this section we generalize our construction to an 
arbitrary even number N — In of subsystems. We use 



the notation cr 



(N) 



1 xp — ®x. . .xp. . .p 

tion of the operator x± . . . x n p n +i . . .p2n and, similarly, 
Opx — cr p ... px ...x stands for the standard deviation of the 
operator p\ . . .p n x n +i ■ ■ ■ x<z n - The following equalities 
take place: 



for the standard devia- 



JN)2 = 
xp 



A N ) 2 - 



R 



d n tp 



™ \dx n +i . . . dx 



N 



dx, 



d n ip 



dxi . . . dx r , 



(119) 



Here we again consider only spherically symmetric wave 
functions ip(xi, . . . ,xn) — vE'(r), where, as before, we de- 



fine r = j|xj| = y/xj + ... + x' 2 N . 
fives we have 



For the partial deriva- 



dxi 



r dr 



and, in general, for distinct indices i%, 



(120) 
, ik we have 



ld_ 

r dr 



(121) 



For a spherically symmetric wave function Eq. ( 45 ) can 
be generalized as follows: 



where the right-hand side is given by the integral 

2 



(ZW) = , x 
'r 



. x 



N 



r dr 



(122) 



*(r) dx. (123) 



To evaluate this integral we the spherical coordinates 
given by 



(124) 



x\ = r cos (/?!, 
X2 = r sin ipi cos ip2 , 
x 3 = r sin ipi cos Lp 2 cos ip^, 

Xn-i = r sin <p\ sin </? 2 . . . sin ifiN-2 cos ifff-i, 
xn = r sin <p\ sin<^2 ■ ■ ■ sintpN-2 sinipN-i, 

where ipi, . . . , <Pn-2 € [0, tt) and <Pn-i G [0, 2tt). The 
Jacobian of this variable transformation is 

j N—l ■ N-2 ■ N—3 ■ 2 

J=r sin (pi sin tp2 ■ ■ ■ sin tpN-3 sin (pN-2- 

(125) 

We then have the following equality: 



cosVsin^-^-V 



cos 2 ^sin 2 ^" 2 ^" 3 ^ 



. cos 2 (/?at_i sin 2 (Pn-i- 



(126) 



The integrals over the angles can be easily taken explic- 
itly pS] 



sin™ ip cos" 1 ipdip = B 



n + 1 m + 1 



(127) 



2 ' 2 

where B(x,y) is the Euler beta-function. We can write 

'3 3(JV-1)\ 

(128) 



x 2 . . .x% J dip = B ( , 



3 3(JV-2) 



,5 



3 3 • 2\ tt 
2' ~2~ 



where f2 = [0, it) 



N-2 



[0, 27r). Using the well-known 



relation for the beta-function 



B(x,y) 



i»r( y ) 

r(i + y) ' 



(129) 



the right-hand side of Eq. ( 128 ) can be simplified as fol- 
lows: 



x\ . . . x^J dip 



rN 



(130) 



This expression is valid for all N ^ 1, including odd 
values. In our case N — 2n is even and this expression 
can be further simplified: 



n X l--- X lnJdip- 2 2„-l (3n _ 1)r 



Finally, for the averaged value (Z^ N '} we have 



(131) 



22n-l(3 n _l)| 

-\-oo 



1 d 

r dr 



(132) 



*(r) dr. 



For n — 1 this expression coincides with Eq. (41 1. From 
the normalization condition 



^> (xi, . . . , xn) dx = 1 



(133) 



R" 



in the same way we can obtain the following equality 
2tt" 



(n 1)! Jo 



„2n-l, T ,2 



* 2 (r) dr = 1. 



(134) 



For n — 1 it coincides with Eq. (42 1. Let us introduce 
the function f(r) via the relation 



*(r) = V^/(r 2 "). (135) 

V 7T" 

which generalizes Eq. (46). From Eq. ( |134 ) we have 

2tt" n! /- +oc 
(n-l)!^7 



2/ 2n\ 7 2n _ / i*2 



(136) 



f z {r 2n )dr 2n = / r(r)dr = 1. 
o 
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We see that the function f(r) is normalized in the ordi- 
nary sense. Computation of the integral ( 132 ) requires 
more work. As one can easily check, we have 

1 d 

) f[T") = 



dr 



f(r 2n ) = J2 a k r2n(k ' 1) f {k) (r 2n ), (137) 



In fact, the elements of these matrices are 



Aij — 



i < j 





(-i)' +J 



The matrix element of the product is 



i < j 

i > i 
(144) 



where the coefficients a k are to be determined. We can 
write the Taylor expansion of the function f(r 2n ) as 



g /(fc)(0) 
n=0 



2nk 



fc! 



(138) 



and apply the differential operator to this expansion 



ld_ 
r dr 



g/W(0) 2 n (nk)! 
fc=i 



fc! (n(fc-l))! 



„2n(fc-l) 



_ 2 ny / (fe+1) (0)4fc + l)! 
A. (fc + i)i {nk)[ ' 



,2nk 



(139) 



On the other hand, each of the n derivatives on the right- 
hand side of Eq. ( 137 ) can also be expanded into its Tay- 



lor series and when these expansions are substituted into 



Eq. |137j) we get 

n +og 



k=l j=o 



3- 



(140) 



The two expressions ( 139 ) and ( 140 1 must agree, so 
the coefficients of the powers r 2nK must be equal for 

all fc = 0, 1, When we equate these coefficients for 

fc = 0, 1, . . . ,n — 1 we get the following system of equa- 
tions: 



E 



2 n 



(nm 



2 n nl 



(141) 



-J (to — fc)! m\(n{m — 1))! ml \ n 
for fji — 1, . . . , n. This system has the following matrix: 



.4 



I * 
v 

2! 





1 

1! 



\(n-l)! (n-2)! (n-3)! 






.. 0\ 





.. 


1 

0! 


.. 


1 

-3)! ■ 


■■ h) 



(142) 



One can easily verify that the inverse matrix A 1 is given 
by a similar expression 



A- 1 = 



2! 





_°i 

1! 













\ (n-l)l (n-2)! (n-3)! 



(143) 



k) 



{A ■ A^ 1 )^ = ^2 A lk (A~ 1 ) k 



kj- 



(145) 



fe=0 



Only the terms with i ^ k and k j are nonzero in this 
sum, so all terms above the main diagonal, i.e. the terms 
with i < j, are zero (the product of two lower triangular 
matrices is also lower triangular). For i j we have 



(A-A-% = J2 



k=3 



1 (-l) fc +J 

(i-fc)! (k-j)\ 



y ■" k=0 



1-3 



(146) 



kl l -3 
k 



= 6 



*3> 



where the last equality is due to the well-known identity 
for the binomial coefficients [17] . We see that the product 
of these two matrices is the identity matrix, A-A^ 1 = E n , 
so the matrix A^ 1 defined by Eq. (143) is indeed the 
inverse of the matrix ( 142 ). 



We now can express the coefficients a k explicitly 

k 

ak 



^-^ (h — ill j\ 



\ (k-jy 
(-i) 



72 



fc! 



k k 

-EM) 



(147) 



Now one can verify that coefficients of the powers 



„2nk 



0,1,. 



in the expression (139) and (140) agree not only for 
k 



1, but for all fc. When we substitute 



Eqs. ( 135 ) , ( 137 1 into Eq. ( 132 ) we get an explicit expres- 



sion for (Z^) in the form of the following functional: 

2 



£(n)\ = 3 



(n!) 



P|3 



(3n)! 

where the coefficients b k read as 

^k 



,.+00 / n \ z 



bk = 



(~1) A 
fc! 



3 = 1 



(148) 



(149) 



For n = 1 Eq. (148) coincides with Eq. (50). 

First , we show that the integral on the right-hand side 
of Eq. ( 148 1 is never zero for a square-integrable function 



f(r). In fact, if this integral is zero then f(r) satisfies 
the following differential equation of the nth order on 
the interval (0, +00): 



£>r fc /«(r)=0. 



(150) 



fc=i 
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This equation has n linearly independent solutions. We 
can try to find them in the form f(r) = r Q , where a must 
satisfy the equation 



J2*>k(a) k = 0, 



(151) 



fe=i 



where (a)k = oc{a — 1) . . . (a — k + 1) is the Pochhammer 
symbol. One can easily verify that this equation has n 
different roots 

1 n- 1 



(152) 

This follows from the fact that if we take f(r) — r^ n , 
then f{r 2n ) = r 2k and 



1 d\ „, 2n\ _ d \ r 2k _ g 



r dr I \rdr 



(153) 



for k = 0, 1, . . . , n — 1. In Appendix [C] we prove directly 
from the expression for the coefficients bk, Eq. ( 149 ), that 



the numbers ( 152 ) satisfy the equation ( 151 1. We see that 



any solution of the equation (150) has the form 



f(r) = C„_i v^i + . .. + Cx^r + Co 



(154) 



None of these functions is normalizable on the interval 
(0,+oo). 

Now we find the infimum of the operators (Z^ n ') over 
the set of all normalized functions f(r). It is rather a 
challenging task to do it in general so we find the min- 
imum for n = 2, 3 and develop a method which allows 
to do it in general. In both cases the functions obtained 
in the previous step are used. In contrast to the bipar- 
tite case, here we could not find explicit analytical ex- 
pressions for the solutions we obtain, but we developed a 
technique to prove that these solution are indeed optimal 
without having explicit expressions for them. 



A. Four-partite case 

For N — 4 we need to minimize the integral 

= - jf (rf(r) + 2r 2 /" (r)) 2 dr. (155) 



We prove a more general 
Theorem 3. The equality 



inf \\rf + ar 2 /'T 2 



1 /3 



-a- 1 



(156) 



||/||=i ■■ - 4 V2 

is valid provided that a ^ 1 . 

Another formulation of this theorem is the statement 
that for an arbitrary non-normalized (but normalizable) 
function / we have the inequality 

\\rf' + ar 2 f"\\>± (\a -l)||/||, (157) 



and the coefficient on the right-hand side is the best pos- 
sible (for a ^ 1). 



Proof. Computing the scalar product 



(r/'W/V/')=(l-^J Wrf'f, 
we can apply Cauchy-Schwarz inequality to get 



\\rf+ar 2 f"\\ > [ ^ - l) \\rf\\ > 



1 

2 V 2 



-a- 1 



(158) 



(159) 



where we used the inequality ( 56 ) and the normalization 
ll/H = 1. We need to show that one can find such func- 
tions / with ll/H = 1 that the difference between the 



left-hand and the right-hand side of the inequality ( 156 ) 



becomes arbitrary small. We construct such functions 
with the help of the functions we obtained for N = 2. 

Let us take the function f%(r) = /^ (r) defined by 
Eq. (97) and find a normalizable solution g^(r) of the 
equation 



(1 - a)g c + arg's = fa 



(160) 



To solve it, we first solve the following auxiliary equation, 
motivated by the integral representation |96|): 



(1 - a)G + arG' 



e" 7r , 7 



1 1 - V^cos6> 

2 1 + v^cos0' 



Then, the corresponding solution of Eq. ( |160[ ) is 
1 r G(r) 



.%0) = 



y/2nK(X) J a 1 + V?cos0 



de, 



(161) 



(162) 



where, in fact, G(r) depends on both £ and 9, since 7 
in Eq. ( 161 1 depends on them. The general solution of 
Eq. qi61b is 



G{r) = cr~ 



1 



~(7r) ■ T 



a - 1 



,7r 



(163) 



so a normalizable solution of Eq. ( 160[ ) is given by 

1 r ( 7 0^r(-^i, 7 r) 



.9?0) = - 



« v / 2^(£) Jo 



1 + a/? cos 



where r(£, x) is the incomplete Gamma-function 



r+00 

r(c,x) = / e-Ht-t-dt. 

J X 



d0 7 
(164) 



(165) 



The solution ( 163 1 for c = can be also written as [TH] 



G(r) 



ar(2 



+°° e~H c - 



■yr + t 



dt, 



(166) 



so it is clear that the corresponding solution ( 164 1 is nor- 



malizable. We could not find an explicit expression for 



14 



the function <?£ (r) itself or for its norm in terms of known 
special functions, but it is possible to find the limit of 
H^ll when f -> 1. 
From Eq. ( 160 1 we obtain 



rg' 6 + ar 2 g'l=rf' s , (167) 
and, since the solution g^(r) is normalizable, we have 

(1 - o)(l - 2a)||^|| 2 + a 2 \\rg>\\ 2 = ||/ e || 2 = 1, 
(l-3a)\\rg^ + a 2 \\r 2 g'^ = \\rf^. 

To estimate ||r 2 <7^||, let us compute the scalar product 

3 



(r g'l,rg'(.) 



\ r 9e,\\ , 



from which we get the inequality 

3 



r 2 9'l\\>~\\r9' £ l 



(169) 



(170) 



and from the second equality of Eq. ( 168 1 we get 
9 



rf'f^[l-3a+ja')\\rg' s r=[^a-l\ \\r£ 



I 1 1 2 



J 1 1 2 



(171) 



Then, from this inequality, the first equality of Eq. ( 168 ) 
and the inequality ( 56 ) we have 



1 < (4(l-a)(l-2a) + a 2 )||r^|| 2 

'3 V (172) 
ja-lj M| 2 ^4||r4|| 2 . 

It is in this place that we use the condition a 1 that 
guarantees non-negativity of the coefficient (1— a)(l— 2a). 
By construction of the function /t, we have 4|jr/^|j 2 — > 1 
when £ — > 1, from which we immediately get 



I 1 1 2 



4(§a-l) 



(173) 



and thus \\gt\\ — > (§ a ~ l) when £ — > 1. Results of 
numerical integration of g^ for several randomly chosen 
a > 1 and for £ 1 agree with this conclusion (within the 
accuracy provided by the numerical integration routines). 
If we take the function 



gd r ) 
Ikll 



then this function is normalized, \\g^\ 
(0, 1), and we have 



(174) 
1 for all £ € 



m+ar 2 ~g'l\\ 2 = 



\\rf'A\ 2 



1 /3 



IflCir 4 V2 
when £ — > 1, which concludes the proof. 



-a - 1 



(175) 
□ 



From this theorem, for a = 2, we get inf(Z( 4 )) = 1/30. 
As a minimizing family of functions we can take the func- 
tions jj^{r) — <?f(V)) where <?g(r) is given by Eqs. fl74| 



and (164 1 for a = 2. The corresponding function 

Eq. jl35| , is shown in Fig. [ij The wave function ip(x) is 

given by 



\/2 (2) , , 2 . 2 i 2 i 2\2\ 
= /f '((Zi + Z 2 + z 3 + X 4 ) ). 

7T s 



(176) 



The functions (174) for a = 3/2, which we will denote 

(4} 

by the same symbol (r), will be used in the next 
subsection to construct a minimizing family of functions 
for the case N — 6. 




FIG. 4: The wave function \&(r), defined by Eqs. (1351 for 



n = 2, ( |174[ ) and ( |164[ ) with a = 2 for a few different values 
off. 



B. Six-partite case 

For iV = 6 we have to minimize the integral 

f+oo 



1 / 

<^ (6) > = ^ / [rf(r) + 9r 2 f"(r) 



560 



n 



,2 (177) 
+ ~r 3 /"'(r)J ^. 

The solution is given by 

Theorem 4. The following equality holds true: 

35 



mf(Z^) 



4096 



(178) 



^4 family of functions that approach this infimum can be 
constructed with the help of the functions obtained for the 
previous case n — 2. 

Proof. We have 



r/' + 9r 2 /"+ V/'" 

= r(3r/'-2/)' + ^r 2 (3r/'-2/)", 



(179) 
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and from Eq. (157) we get the inequality 



\rf> + ^f» + \r*n\>\(\-\ - 1)11-/ -*'f'\ 



= ^||2/-3r/'|| 



(180) 



The value ||2/ — 3r/'|| 2 can be estimated as follows: 



|2/-3r/|| 2 = 10||/|| 2 + 9||r/' 



> 10 



9 



\f\\ 



49, 



/II 



49 



(181) 



where we used the inequality ( 56 1 . Combining this result 



with the inequality ( 180 1 we get 



\\rf + 9r 2 f"+ 9 -r 3 n\>l-l 



35 
16' 



(182) 



We now prove that this estimation is the best possible. 
Consider a normalizable solution of the equation 



2hc + '3rti ( = ,/ e l 



(4) 



(183) 



where /^(r) is the function given by Eqs. (174) and 



(164) for a — 3/2. This function fc 4 \r) has the following 



integral representation: 



« ^r(-|, 7 r) 

1 + V? cos 9 



d6, 



(184) 



W II - 



where -/V(£) is the normalization such that ||/^ 
To solve Eq. ( |183[ ) we first find a solution of the equation 



(185) 



2H ( .+ZrH' ( : = ^rFT[--,jr 



and then construct from it the corresponding solution of 
Eq. Jl83} via 



H 6 (r) 



d0. 



(186) 



A normalizable solution of Eq. ( 185 ) reads as 



Hdr) = 2 e 



(187) 



or, according to [19] . this solution can also be written as 



7r 



3 "/re 



-yr f+OO 



2 r(|) Jo Vthr + t) 

- lr ,+oo y~ te ~t 

dt, 



dt 



(188) 



r(|)7 ir + t 



from which it is clear that this solution is normalizable. 

We could not find an explicit expression for the func- 
tion h^(r) itself or for its norm ||/if|| but we can prove 
some statements about the behavior of the norm II he II 



when £ — > 1. From Eq. (183) we have 

rh't + 9r 2 /^' + 9 -r 3 hf = rff + \r 2 ff" . 



(189) 



Taking the norm of the both sides of this equality, we get 



28|K|| 2 -^| 



r 2 /i'"' 2 



— I 

4 1 



r 3 h'f'\\ 2 



(4)/ 



3 r 2 ,(4)// 



(190) 



where we used the relations 



(rh' v r 2 h'i) = -^\\rh' £ \\ 2 , (rX',^") = - 
(rh' s ,r s hf) = 6\\rh' s \\ 2 -\\r 2 ti{\\ 2 . 



r 2 h"\\\ 



(191) 



From the second of these relations we can get the follow- 
ing inequality: 



(192) 



If we substitute it into the left-hand side of Eq. ( 190 ) we 

get 



28||r/ l '|| 2 -^||r 2 / l ''|| 2 , 



63 
16 1 



(193) 



so the second coefficient is negative and we cannot fol- 
low the idea of the prove for n = 2. Up to now we 
used relations like ( 192 ) which are derived from Cauchy- 



Schwarz inequality and this inequality is based on the 
non- negativity of the expression \\af + f3g\\ 2 for all real 
numbers a and /3. Here we need to go one step fur- 
ther and use the non-negativity of similar expression with 
three terms. Namely, we can write 



arh' e + (3r 2 ti i ! + -r 3 ti e " 



(a 2 -3a/3 + 54a) || rk 



I ||2 



81, 



3l///||2 



^ 0. 



Then we find such a and f3 that 

a 2 - 3a/3 + 54a = 28-49 



f3 2 -9a- = - 



261 



(194) 



(195) 



2 ' 2 

One can verify that this system has two solutions 

a=^(9±V74), 0= ^(24 ±y/U). (196) 
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If we take one these two solutions, (ao,/3o), the equality 
( 190 1 can be written as 



49 



\rh'\\ 2 



a rh^ + p r + -r 



Ai)> i ^ 2 A<t)r> 



so that we have the inequality 

2 



49 



2 Ai)" 



From the differential equation ( 183 ) we have 
10||/i e || 2 + 9||r/4|| 2 H|4 2) || 2 = l 



for all £. Then, using Eq. (198), we have 



1 = W\\h ( \\ 2 +9||r^|| 2 49||r/4| 



Sl 5 



A 4 )' i " * t 



2 /(4)» 



We have chosen the function fj^ (r) in such a way that 



(197) 



(198) 



(199) 



(200) 



The corresponding function \&(r) is shown in Fig. [5] and 
the corresponding wave functions read as 



ip(xi, ...,x 6 ) = #(r) 



(206) 

As in the previous case, we have not been able to find an 
explicit expression for /^ 6 ' (r) in terms of known special 
function, but we could prove that these functions have 
the desired property without such an expression. 




FIG. 5: The wave function *(r), defined by Eqs. (135} for 
n = 3, (|202|, (fl86l and 11871 for a few different values of £. 



2 M)» 



1 /3 3 
4 \ 2 ' 2 



(201) 



when £ 1. From Eq. (200) we get that ||r/il|| -> 1/7 



and then from Eq. (199) we obtain that 1 1 \ \ — > 2/7 when 
£ — > 1. Results of numerical integration for £ w 1 agree 
with this conclusion. Now, let us take the function 



fe c (r) 



(202) 



This function is normalized, \\h^\\ — 1, and from 
Eq. (p9| we get 



r% + 9r 2 /4' + r 3 ^" 



1 




INI 


7 5 


35 


-> - • - 




2 8 


~ 16 



"V f (4)// 

2 ; « 



(203) 



when £ — > 1. From the definition of Z^ 3 \ Eq. (177), we 
have 



35 
4096' 



which concludes the proof. 



(204) 
□ 



The functions f^ e \r) that minimize (Z^) can now be 
taken as 



jf(r) = ^(r). 



(205) 



VII. CONCLUSION 

In conclusion, we summarize our results obtained in 
this work. We have tried to minimize the uncertainty 
product a xp a px and its multipartite generalizations. We 
have shown that pure states are sufficient to minimize 
such uncertainty products. We have outlined the gen- 
eral approach for the case of an arbitrary even number of 
parties and explicitly found the infimum of these prod- 
ucts over a special class of states with real spherically 
symmetric wave functions in bipartite, four-partite and 
six-partite cases. We have also constructed parametrized 
families of states that approach the corresponding infi- 
mums by varying the parameter. Namely, we have shown 
that any bipartite state for which 



1 



1 



(207) 



is entangled; that any four-partite state for which 

^ @xxpp®ppxx ^ (208) 
is entangled; and that any six-partite state such that 

(209) 



35 1 _ 35 1 
64 ' 64 ~ 4096 < a xxx ppp° ppp xxx < ~ 



64 



is entangled and provided examples of states that tend 
to the limits of 1/8, 1/30 and 35/4096 respectively. The 
next step will be to answer the question whether these 
limiting values can be improved by using more general 
states. 
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Appendix A: Changing the order of integration 

We need to verify that we can exchange summation 
and integration in Eq. ( 77 ) . We use the following theorem 



that can be found, for example, in [20] . 

Theorem 5. Assume that the sequence of functions 
F n (t), n = 0,1,..., defined on an infinite interval 
[tQ,+oo), uniformly converges to the function F(t) on 
any finite interval [to,T], T > tQ. If the integral 
J t ^°° F n (t) dt converges uniformly with respect to n, then 
the integral J t + °° F(t) dt exists and the equality 



lim 



F n (t)dt= / F(t)dt. (Al) 

/to « to 

is valid. In other words, under these conditions the limit 
of integrals is equal to the integral of the point-wise lim- 
iting function. 

We apply this theorem to the functions 
c e r S n (tt)J (2Vr~t)e- t 

c e Vte b n {€,t) , 



F n {t) 



(A2) 



where £ and r are fixed. From the representation ( 82 ) it is 



clear that the partial sum S n (£, t) converges uniformly to 
Iq (iy 7 ?) on any finite interval. To prove that the integral 
Jo + °° CO conver g es uniformly note that 

(X^te-'Sni^t) < y/te-*I Q {ty/l). (A3) 

The Bessel function io(0 has the following asymptotic 
expansion for large t |21) : 

MO ~ -£= ( 1 + 4 + + • • )> (A4) 



from which we immediately get that v / ie~*/o(iv / C) 
when t —> +oo (take into account that < £ < 1) and, 
as a consequence, we see that the function v / ie~ t / (i-\/£) 
is bounded on [0, +oo). We have just proved that the 
functions \/te~ t S n (^,t) are uniformly bounded. On the 
other hand, the function Jo(2y/rt) / y/i (which does not 
depend on n) is integrable for r > [22 



Vt 



dt = 2 



J (2ty/f)dt= (A5) 



From this we can conclude that the integral J +O ° F n (t) dt 
converges uniformly and thus we can apply the Theo- 
rem 1 to verify the correctness of the expression ( 83 1 . 



Appendix B: Normalization of the coefficients c nm 

Here we prove that the normalization condition for the 



coefficients defined by Eq. ( 108 1 

+oo 



E 

n,m— 



c z =1 



(Bl) 



can be obtained directly from the explicit expression 



(118) for these coefficients. In fact, we have 



E 



c 2 = 



E E < 

N=0n+m=N 



E' C nm, (B2) 



N=0n+m=AN 



where the prime means that the sum runs over only those 
pairs of n and m that have the form n = An' , m = Am' 
or n = An' + 2, m = Am' + 2. The first equality is simply 
the " diagonal" summation of this double series with non- 
negative coefficients (so the sum does not depend on the 
order of summation) and the second equality is valid since 
if N is not a multiple of 4 then c nm = for all pairs (n, m) 
with n + m = N. From the expression (118) for the 



coefficients we see that we need to compute the following 
primed sum: 



S N = E ( " 

n+m=4N * 2 



N 

E 

71=0 



An\ ( AN - An 



2nJ \2N -2n 



N-l 

E 



(B3) 



An + 2\ /AN -An -2 
2n + 1 1 [ 2N - 2n - 1 



To compute this sum, let us introduce the function F(x) 
via the equality |23j 



+00 



F&) = E 



2n 



1 



Vl^Ax 



(B4) 



Then we can observe that the sum (B3| can be written 



in a compact form with the help of F(x) as 



27V 



S N = ^[x™]F(a;)[x 2Ar ^"]F(2;) = [x 2N ]F 2 (x), (B5) 
n=0 

where \x n ]F(x) is the coefficient of x n in the Taylor ex- 
pansion of F(x). Thus, we have 



Sn — 



„2/Vl 



T -Ax 



= 2 



4N 



(B6) 



Now we can finish the computation of the sum on the 
left-hand side of Eq. ([Blj): 



E 



+00 /„ Ar \ 2 

2K(£) ^ V N 

vsy N—0 v 



2K(£) ^ 



(2N 

+00 /ri ArX 2 



N 



f2N 

[n 



16 2 

2\ N 



N 



(B7) 



= 1, 



where we have taken into account the relation (72 ). This 



completes the proof of the equality (Bl ) 



Appendix C: Roots of Eq. (151 1 



Here we prove that the numbers ( 152 ) are the roots of 



the equation ( 151 1 where the coefficients bk are given by 
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Eq. (149). We have 



for a given by Eq. (152 1. For a — this is clearly true 



since (0)„+i = 0. For a = i/n, i = 1, ...,n— 1, we have 



= EE 

fe=ii=i 



fe=i 



(-] ^ 



fc! 



(")fe 



E(-d j 



(Cl) 



5 



n I m \n 



n+l j — i 



n-1 



i=0 



where we have exchanged the summationorder accord- gince the binomial coefficien t (J«) can be expanded as 
ing to the equality Y%=i E,*=i = E?=i Efc=,- The inner 



sum on the right-hand side of Eq. (Cl ) is easier to com- 



pute when a is a sufficiently large integer number. In 
this case we can write {ot)k = oil /{a — k)l and transform 
this sum as follows: 



(-i) r 



a — j — 1 
n-j 



(-!)"+! (a) n+1 fn 
nl j -a \j 



(C2) 



where we have used a simple relation for the binomial 
coefficients 1241. Note that 



(a)„+i a(a-l)...(a-n) 



3 



J 



(C3) 



is a polynomial of a for j = 1, . . . , n. Since the polyno- 
mial on the left-hand side of Eq. ( |C2[ ) equals the poly- 
nomial on the right-hand side for all sufficiently large 
integer a, these two polynomials must be same, and the 
equality (C2| must be valid for all a. When we substi- 



tute the expression on the right-hand side of Eq. (C2) 
into Eq. (Cl ) we see that we need to prove that 



5(a) = (a) n+1 ^2 



try 



(C4) 



jn 

n ) n 



n— 1 



(C6) 



and l/(j — — ) = n/(nj — i) cancels one term of this 
product. Since i ^ the term with I — 0, i.e. jn, is 
always present. Then we can expand the product on the 
right-hand side of Eq. ( C5 ) and get 



*(s)-S>S>' 



)^(.j^ = 0, (C7) 

where A p are some numbers and each inner sum is zero 
due to another standard relation for binomial coefficients 



4=1 4=0 VJ/ 



(C8) 



for p = 1, . . . , n — 1. It is important to note here that if 
j ^ then p in the sum ( |C7|) sta rts from 1 and in this 
case the first equality in Eq. (|C8| holds true. 
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